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HIGHLIGHTS 


►  We  developed  a  theoretically  based  thermal  model  for  lithium  ion  battery  pack. 

►  Increasing  current  rate  has  more  obvious  effect  on  irreversible  heat  generation  rate. 

►  The  change  of  SOC  mainly  affects  the  reversible  heat  generation  rate. 

►  There  always  exists  an  optimal  resistance  coefficient  for  battery  pack. 

►  Reversible  heat  affects  battery  temperature  only  during  continuous  charge  and  discharge. 
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Based  on  the  porous  electrode  and  concentrated  solution  theory,  a  thermal  model  is  developed  for 
lithium  ion  battery  pack.  The  accuracy  of  predicted  battery  temperatures  is  validated  by  charge 
—discharge  cycling  experiments  under  natural  and  forced  convection  conditions.  The  heat  generation 
and  dissipation  rates  of  battery  under  different  conditions  are  simulated  by  the  proposed  model  and  the 
results  indicate  that:  (1 )  the  SOC  change  has  a  significant  effect  on  the  reversible  heat  generation  rate  but 
has  almost  no  influence  on  the  irreversible  heat  generation  rate;  (2)  the  generation  rates  of  reversible 
and  irreversible  heat  during  charge  are  almost  equal  to  that  during  discharge  with  the  same  SOC  and 
current  rate,  but  the  effect  of  reversible  heat  on  battery  temperature  is  opposite;  (3)  for  enhancing  heat 
dissipation  with  a  given  input  power  of  cooling  fan,  there  always  exists  an  optimum  value  for  the 
resistance  coefficient  of  battery  pack,  and  the  optimal  coefficient  is  increased  when  the  input  power  of 
fan  increases.  In  addition,  the  comparisons  between  the  predicted  and  measured  battery  temperature 
indicate  that,  the  reversible  heat  has  significant  influence  on  battery  temperature  during  continuous 
charge  and  discharge,  especially  under  low  current  rate,  but  the  influence  can’t  be  observed  during 
charge-discharge  cycles. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  ion  battery  is  considered  as  one  of  the  most  ideal 
energy  storage  equipment  for  electric  vehicles  by  virtue  of  its 
advantages,  such  as  high  power  density,  long  cycle  life  and  low 
self-discharge  rate  etc.  However,  the  temperature  of  lithium  ion 
battery  is  gradually  increased  during  the  charging  or  discharging 
process.  If  battery  temperature  exceeds  the  maximum  allowed 
value,  the  degradation  rate  of  battery  life  would  be  accelerated  and 
thermal  runaway  would  even  be  triggered  under  some  extreme 
conditions  [1].  Therefore,  a  thermal  management  system  which 
can  maintain  battery  temperature  below  the  limiting  value  is 
required  to  be  specially  designed  for  the  lithium  ion  battery  pack 
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used  in  electric  vehicles.  In  order  to  ensure  that  the  battery 
temperature  can  be  effectively  controlled  below  the  limiting  value 
under  any  driving  conditions,  it  is  necessary  to  fully  understand 
the  characteristics  of  heat  generation  and  heat  dissipation  of 
lithium  ion  battery  pack.  Numerical  modeling  is  a  better  way  to 
study  the  thermal  behavior  of  battery  pack  than  experimental  tests 
with  respect  to  the  cost  and  time  consumption.  Therefore,  devel¬ 
oping  a  thermal  model  to  predict  the  heat  generation  rate  (HGR) 
and  heat  dissipation  rate  of  the  single  batteries  within  lithium  ion 
battery  pack  is  helpful  and  crucial  for  properly  designing  the 
thermal  management  system. 

Several  thermal  models  of  lithium  ion  battery  pack  have  been 
presented  in  published  literature.  According  to  the  difference  in 
calculating  the  HGR  of  lithium  ion  battery,  the  thermal  models  can 
be  classified  into  empirically  based,  semi-empirically  based  and 
theoretically  based  thermal  model.  In  the  empirically  based 
thermal  model,  the  battery  HGR  is  determined  purely  by  the 
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experimental  data  of  the  battery  modeled,  e.g.  the  model  devel¬ 
oped  by  Park  et  al.  [2]  and  Mi  et  al.  [3].  Park’s  model  calculated  the 
battery  HGR  by  Joule’s  law  with  the  internal  resistances  determined 
from  fitting  experimental  data  of  battery  at  different  temperatures, 
and  Mi’s  model  deduced  the  HGR  by  experimentally  measuring  the 
temperature  evolution  of  battery.  The  battery  HGR  in  semi- 
empirically  based  thermal  model  is  calculated  by  Bernardi  equa¬ 
tion  [4],  e.g.  the  model  established  by  Chen  et  al.  [5]  and  Zhang  et  al. 
[6].  As  the  Bernardi  equation  was  derived  from  the  processes  of  the 
electrochemical  reaction  and  charge  transfer  occurring  inside 
lithium  ion  battery  [4],  the  semi-empirically  based  thermal  model 
is  different  from  the  empirically  based  thermal  model,  though  the 
parameters  involved  in  Bernardi  equation  are  also  determined  from 
experimental  data.  In  the  theoretically  based  thermal  model  the 
battery  HGR  is  obtained  from  electrode  reaction  mechanism,  e.g. 
the  model  proposed  by  Smith  et  al.  [7]  calculated  the  HGR  of 
lithium  ion  battery  based  on  the  porous  electrode  and  concentrated 
solution  theory. 

Since  the  reliability  of  battery  HGR  is  completely  dependent  on 
experimental  data,  the  application  of  empirically  based  thermal 
model  is  limited  to  specific  battery  packs.  The  semi-empirically 
based  thermal  model  has  wider  application  scope,  but  the  model 
accuracy  needs  to  be  improved  because  of  the  assumption  of 
neglecting  the  spatial  variations  of  parameters  such  as  reaction 
current  density,  active  material  concentrations  etc.  in  electrodes, 
which  could  lead  to  significant  errors  in  estimating  battery  HGR  [8]. 
The  theoretically  based  thermal  model  determines  battery  HGR  on 
the  basis  of  the  universal  electrode  reaction  mechanism  and  the 
consideration  of  the  spatial  variation  of  parameters  in  electrode. 
Therefore,  it  has  the  widest  application  scope  and  the  highest 
accuracy  among  the  three  types  of  thermal  models.  However,  most 
existing  research  is  focused  on  the  theoretical  model  of  single 
lithium  ion  battery,  and  little  work  has  been  done  on  the  devel¬ 
opment  of  theoretically  based  thermal  model  for  lithium  ion 
battery  pack. 

This  paper  aims  to  develop  a  theoretically  based  thermal 
model,  which  can  accurately  predict  the  heat  generation,  heat 
dissipation  and  temperature  rise  of  the  single  cells  within  lithium 
ion  battery  pack.  In  addition,  as  the  reversible  heat  neglected  in 
Smith’s  model  [7]  is  indicated  to  account  for  a  significant 
proportion  of  the  total  heat  generated  inside  lithium  ion  battery 
[9—11],  the  battery  temperatures  with  and  without  considering 
reversible  heat  as  heat  source  are  estimated  by  the  model  and 
compared  to  the  measured  temperatures,  in  an  effort  to  investigate 
the  effect  of  reversible  heat  on  the  model  accuracy  under  different 
operating  conditions. 


2.  Model  development 

2.1.  Lithium  ion  battery  pack 

Fig.  1  shows  the  components  and  arrangement  of  a  typical 
lithium  ion  battery  pack  for  electric  vehicle  application.  The  battery 
pack  mainly  consists  of  N  cells  connected  in  series,  (N  +  1 )  cooling 
channels  arranged  between  the  cells,  deflector  plates  and  cooling 
fan  etc.  Under  forced  convection  condition,  a  negative  pressure 
environment  is  created  inside  battery  pack  by  the  cooling  fan.  As 
a  result,  cooling  air  would  be  drawn  into  the  battery  pack  through 
the  inlet  and  guided  into  the  cooling  channels  by  the  lower  air 
deflector.  With  the  guidance  of  the  upper  air  deflector,  the  heated 
air  would  be  expelled  from  battery  pack  through  the  outlet.  An 
ideal  design  of  the  deflectors  can  evenly  distribute  the  airflow  to 
provide  uniform  pressure  difference  between  the  inlet  and  outlet  of 
each  channel. 

The  thermal  model  is  developed  to  assist  in  designing  the 
thermal  management  system  of  lithium  ion  battery  pack.  For 
a  prismatic  lithium  ion  battery  shown  in  Fig.  1,  a  small  temperature 
gradient  in  battery  can  drive  a  significant  amount  of  heat  to  be 
transferred  out,  because  of  the  big  transfer  area  and  small  transfer 
distance  in  the  thickness  direction  of  battery.  Therefore,  under 
normal  operating  condition,  the  temperature  gradient  inside 
battery  is  generally  small  enough  to  neglect  its  impact  on  thermal 
management.  To  increase  the  calculation  speed  of  model,  the 
following  assumptions  are  made: 

(1)  The  heat  generation  and  temperature  inside  a  single  cell  are 
assumed  to  be  uniformly  distributed,  and  the  thermal  capacity 
and  mass  of  the  cells  are  also  assumed  to  have  no  difference. 

(2)  The  pressure  difference  between  the  inlet  and  outlet  of  every 
channel  is  supposed  to  be  the  same,  and  the  velocity  gradient  of 
the  main  airflow  along  cooling  channel  can  be  set  at  zero. 

(3)  The  inlet  air  temperature  of  every  channel  is  equal  to  the 
ambient  temperature. 

According  to  the  law  of  energy  conservation,  the  temperature 
change  of  cell  n  within  the  battery  pack  shown  in  Fig.  1  can  be 
calculated  by  the  following  equation: 

Cbat^bat  =  Qgen(n)  -  Qdis(n)  -  Qcon(n)  (1) 

where  Qgen(n)  is  the  heat  generation  rate  of  cell  n,  Qdis(n)  is  the  heat 
dissipation  rate  of  cell  n,  Qc0  n(n)  is  the  heat  transfer  rate  between  cell 


Fig.  1.  Schematic  view  of  lithium  ion  battery  pack. 
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n  and  adjacent  cells,  Qcon(n)  =  [2Tbat(n)  -  Tbat(n  +  1)  -  Tbat(n  -  1)]/ 
Rcon,  Rcon  is  the  thermal  resistance  of  the  material  located  between 
two  cells.  For  a  given  lithium  ion  battery  pack,  the  temperature 
change  of  cell  n  mainly  depends  on  Qgen(rc)  and  Qdis(n)  as  expressed 
by  Eq.  (1 ).  In  the  following  sections,  we  will  detail  the  mathematical 
expressions  for  them. 


Due  to  the  high  electronic  conductivity  of  the  current  collectors 
bound  to  electrodes,  the  lithium  ions  and  electrons  are  always 
transported  through  the  path  which  is  perpendicular  to  the 
collectors  [12].  Therefore,  we  only  consider  the  spatial  variation  of 
parameters  along  the  thickness  direction  of  electrodes  in  this  paper. 
Without  loss  of  generality,  the  following  assumptions  are  also  made 
for  the  single  cell: 


2.2.  Heat  generation  within  a  single  battery 

The  main  components  of  a  lithium  ion  battery  are  shown  in 
Fig.  2,  which  include  negative  current  collector,  positive  current 
collector,  negative  electrode,  positive  electrode  and  an  electron 
blocking  separator.  Generally,  the  negative  current  collector  and 
positive  current  collector  are  made  of  cuprum  and  aluminum 
respectively,  and  the  electrodes  and  separator  are  made  of  porous 
material  and  filled  with  electrolyte  composed  of  organic  solvent 
(ethylene  carbonate/dimethyl  carbonate)  and  lithium  salt  (LiPFg). 
As  can  be  seen  from  Fig.  2,  LixC6  and  LiMn204  are  separately  used 
as  active  material  for  the  negative  electrode  and  positive  elec¬ 
trode,  and  <5n,  <5P ,  <5sep  are  used  to  represent  the  thickness  of 
negative  electrode,  positive  electrode  and  separator  respectively. 
According  to  the  porous  electrode  theory,  electrodes  can  be 
divided  into  solid  phase  and  electrolyte  phase.  During  discharge, 
lithium  ions  are  extracted  from  the  surface  of  solid  LixC6  particle 
and  electrons  are  released  at  the  same  time,  then  the  lithium  ions 
travel  via  diffusion  and  migration  through  separator  to  the 
positive  electrode,  and  the  electrons  collected  at  negative  current 
collector  flow  into  external  circuit  to  form  discharge  current.  At 
last,  the  lithium  ions  are  intercalated  into  the  solid  LiMn204 
particles  with  the  electrons  delivered  from  the  external  circuit. 
The  transport  processes  of  lithium  ions  and  electrons  during 
charge  are  inverse  to  that  occurred  during  discharge.  Lithium 
ions  extracted  from  the  surface  of  LiMn204  particles  intercalate 
into  LixC6  particles  with  the  electrons  flowing  into  negative 
electrode. 

The  spatial  variation  of  parameters  such  as  reaction  current 
density,  active  material  concentration  etc.  in  electrodes  is  mainly 
caused  by  the  transport  processes  of  lithium  ions  and  electrons. 


sn 


yscp 


Negative 
#'  electrode 

vV.V 


O 

73 

cd 

Q. 

<D 

tn 


••••• 

•••! 

Discharge  » -^*4; 


Charge 


X. 


g  Positive  y 
electrode  I 

••••••I 


A  \  \  / 

LirC6  Electrolyte  LiMn204 


Fig.  2.  Schematic  of  lithium  ion  battery. 


(1)  Active  material  is  composed  by  spherical  particles  with 
equal  radius  and  the  particles  are  uniformly  distributed  in 
electrodes. 

(2)  The  diffusion  coefficient  of  lithium  ion  in  solid  particles  doesn’t 
change  with  the  lithium  ion  concentration  and  can  be  viewed 
as  a  constant  coefficient. 

(3)  The  changes  of  the  volume  and  porosity  of  the  electrodes 
during  charge  and  discharge  can  be  ignored. 

(4)  Side  reactions  inside  battery  can  also  be  neglected. 

Previous  research  results  indicate  that  the  heat  generated  inside 
lithium  ion  cell  can  be  classified  into  reversible  heat  (also  named  as 
reaction  heat),  ohmic  heat  and  polarization  heat  [11,13-16].  By 
using  the  bottom  point  of  the  interface  between  negative  electrode 
and  current  collector  as  the  origin,  and  the  thickness  direction  of 
electrode  as  the  positive  direction  of  x  axis,  a  one  dimensional 
coordinate  was  established  as  shown  in  Fig.  2.  Based  on  the  pre¬ 
sented  assumptions,  the  computing  method  for  the  generation  rate 
of  the  three  types  of  heat  is  discussed  in  the  following  subsections. 
It  is  worth  to  notice  that  heat  generation  rate  with  positive  sign 
indicates  heat  release,  and  that  with  negative  sign  indicates  heat 
absorption. 

2.2.1  Reaction  heat 

Due  to  the  difference  between  the  Gibbs  free  energy  of 
reactant  and  product,  a  part  of  energy  is  always  needed  to  be 
released  or  absorbed  in  the  form  of  heat  to  maintain  the  energy 
balance  of  the  whole  reaction.  This  part  of  heat  is  called  reaction 
heat  and  its  generation  rate  in  the  cell  shown  in  Fig.  2  can  be 
calculated  by 

i 

Qr  =  -Se|jTbatAS/(neF)dx  (2) 

0 

where  j  is  the  bulk  density  of  reaction  current  in  electrode,  j  >  0 
means  lithium  deintercalation  while  j  <  0  means  lithium  interca¬ 
lation.  AS  is  the  entropy  change  of  electrode  material  for  a  reduc¬ 
tion  reaction,  and  can  be  obtained  from  the  equation  AS  =  (all/ 
0Tbat)nF,  where  U  is  the  equilibrium  electrode  potential,  which  is 
a  strong  function  of  the  SOC  and  temperature  of  battery.  The  rela¬ 
tionship  between  U  and  SOC  can  be  determined  by  fitting  experi¬ 
mental  data  at  a  certain  temperature,  while  the  value  of  0U/0Tbat  is 
often  determined  by  measuring  the  variations  of  U  with  battery 
temperature  for  a  given  SOC  [13].  Fig.  3  provides  the  values  of  dU/ 
6Tbat  when  lithium  intercalation  reaction  was  carried  out  at 
different  stoichiometries  6  in  positive  electrode  (with  LiMn204 
active  material)  and  negative  electrodes  (with  LixC6  active  material) 
[11  ].  As  can  be  seen  from  Fig.  3,  0U/0Tbat  changes  significantly  with 
6.  6  is  in  one-to-one  correspondence  with  the  battery  SOC,  which 
can  be  expressed  as  SOC  =  100%  x  {8  -  #o)/(#i  -  #o),  where  8\  is  the 
stoichiometry  at  100%  SOC  and  do  is  the  stoichiometry  at  0%  SOC. 
The  value  of  6  is  equal  to  cse/csmax,  where  cse  is  the  lithium 
ion  concentration  on  the  surface  of  the  solid  particles.  Using 
the  particle’s  center  as  origin  to  establish  a  spherical  coordinate, 
the  dynamic  change  of  cse  during  charge  and  discharge  can  be 
written  as 
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Fig.  3.  Variations  of  dUldThat  with  stoichiometry. 
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(3) 


with  boundary  conditions  (Ocs/0r)|r=o  =  0,  (0cs/0r)|r=J?  =  -j/ 
(DsFSe). 


2.2.2.  Ohmic  heat 

Electrons  and  lithium  ions  would  collide  with  the  transport 
media  in  the  solid  and  electrolyte  phase  during  their  transport 
processes.  The  heat  generated  by  the  collisions  is  called  ohmic  heat 
and  its  generation  rate  is  proportional  to  the  transport  fluxes  of 
electron  and  lithium  ion.  As  described  by  Ohm’s  law,  the  charge 
transport  in  solid  phase  can  be  expressed  by 


d(Pe 

dx 


dx  +  I2Rc 


(7) 

2.2.3.  Polarization  heat 

Under  open  circuit  condition,  the  activation  energy  barrier 
needed  to  be  overcome  for  deintercalation  reaction  is  equal  to  the 
one  needed  for  intercalation  reaction.  Only  when  the  balance  of 
activation  energy  barrier  is  broken,  can  the  net  current  be  produced 
at  the  interface  between  solid  phase  and  electrolyte  phase.  The 
energy  used  to  break  the  barrier  balance  during  charge  and 
discharge  would  be  released  in  the  form  of  heat  at  last.  The  heat  is 
called  polarization  heat  and  can  be  expressed  by 


l  i 

Qp  =  SeJ  jy dx  =  Se  J j(<ps  -<pe-  U) dx  (8) 

0  0 


where  p  is  the  activation  overpotential.  The  change  of  activation 
energy  barrier  can  be  described  by  rj.  For  example,  rj  >  0  indicates 
the  barrier  needed  to  be  overcome  for  deintercalation  reaction 
would  decrease  aaneFri  and  the  one  needed  to  be  overcome  for 
intercalation  reaction  would  increase  acneFr],  while  r\  <  0  means  the 
former  one  would  increase  -aaneFri  and  the  later  one  would 
decrease  -acneFrj.  rj  =  0  indicates  the  barriers  are  unchanged.  aa 
and  ac  are  the  anodic  and  cathodic  transfer  coefficients  respec¬ 
tively.  The  relationship  between  rj  and  j  can  be  written  as 


j 


seio 


exp 


^bat 


KseiA 
Se  J 


-  exp 


-OcF 

R^bat 


KselA 

Se  J 


(9) 


where  io  is  the  exchange  current  density  and  can  be  determined  by 

l’o  =  /<F(Ce)aa(Cse)Q!c(Csmax  —  Cse)^3  (10) 


j  = 


A  /AeffdgA 
Ox  V  dx  J 


(4) 


with  boundary  conditions  —  creff  (6<^s/ax)  |x=0  =  oe{{(d(ps/dx)  \ 
x=i  =  //Se,  (d(ps/dx)\x=6n  =  {dcps/dx)\x=L_8p  =  0,  where  I  is  the 
working  current  of  battery.  /  >  0  indicates  charge  while  I  <  0 
indicates  discharge.  Due  to  the  lithium  ion  diffusion  caused  by  the 
concentration  gradient  in  electrolyte  phase,  Ohm’s  law  cannot  be 
directly  used  to  calculate  the  transport  flux  of  lithium  ion.  It  can  be 
obtained  by  the  following  expression  derived  from  concentration 
theory  [17]: 


_0^ 

Ox 


+J  =  o 


(5) 


with  boundary  conditions  (0<pe/0x)|x=o  =  (O<Pe/0x)lx=L  =  0,  where 
ce  is  the  lithium  ion  concentration  in  electrolyte  phase.  The  dynamic 
change  of  ce  during  charge  and  discharge  can  be  expressed  as 


0(eeCe)  _  3  fneffSCeA  1  “  f?.- 
at  "  ax  V  e  dx)  +  f  ]  w 

with  boundary  conditions  (0ce/0x)|x=o  =  (0ce/6x)  \X=L  =  0. 
Besides,  the  heat  generation  due  to  the  contact  resistance  Rc 
between  current  collector  and  electrode  can  be  written  as  I2RC 
[7,18].  Therefore,  the  total  generation  rate  of  ohmic  heat  inside 
lithium  ion  battery  can  be  computed  by 


As  discussed  in  subsection  2.2.1,  there  is  a  one-to-one  corre¬ 
spondence  between  the  cse/csmax  and  battery  SOC.  Thus,  the 
exchange  current  density  at  the  interface  of  the  solid  phase  and 
electrolyte  phase  is  closely  related  to  SOC  according  to  Eq.  (10). 

2.3.  Heat  dissipation  of  a  single  battery 

Cooling  air  is  drawn  into  lithium  ion  battery  pack  to  cool  the 
single  cells  through  the  suction  effect  created  by  cooling  fan.  When 
the  air  flows  past  cooling  channel,  a  thermal  boundary  layer  (TBL) 
would  be  formed  at  the  channel  wall  as  shown  in  Fig.  4.  Air 
temperature  in  the  TBL  changes  rapidly  from  the  wall  temperature 
(cell  temperature)  to  the  temperature  of  main  airflow,  and  the 
distance  from  the  wall  to  the  point  where  the  air  temperature  is 
99%  of  the  main  airflow  temperature  is  defined  as  the  thickness  of 
TBL.  Due  to  the  temperature  gradient  established  in  TBL,  heat 
transfer  will  occur  between  the  channel  wall  and  the  air  adjacent  to 
it.  According  to  the  variations  of  the  thickness  of  TBL  and  the  heat 
transfer  in  TBL,  the  cooling  channel  can  be  divided  into  thermal 
entrance  region  and  thermally  development  region  [19].  In  the 
thermal  entrance  region,  the  TBL  is  growing  thicker  and  the  local 
heat  transfer  coefficient  is  decreased  along  channel,  while  in  the 
thermally  developed  region  the  TBL  occupies  the  entire  flow  area 
and  the  local  heat  transfer  coefficient  becomes  constant.  Consid¬ 
ering  that  the  geometric  feature  of  the  rectangular  cooling  channel 
is  similar  with  that  of  a  flat  plate,  the  temperature  profile  in  TBL  in 
the  thermal  entrance  region  can  be  approximated  by  a  cubic 
polynomial  [19].  With  the  2D  coordinate  shown  in  Fig.  4,  the  air 
temperature  in  TBL  can  be  described  by 
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Thermal  entrance  region  j  Thermally  developed  region 


That  («) 

Fig.  4.  Development  of  the  thermal  boundary  layer  in  cooling  channel. 


TayCy)  =  ay3  +  by2 +cy  +  d  (11) 

with  boundary  conditions  Tiay(0)  —  Tbat,  (92Tlay/ay2)  |  q  =  0, 
Tiav(<5)  =  Ta ir,  (0Tlav/0y)|  =  0.  Substituting  Eq.  (11)  into  the  four 

y  I  y=5 

boundary  conditions,  we  get  a  =  (Tt,at  -  Tair)/(2<53),  b  =  0, 
c  =  3(Tair  -  Tbat)/(2<5),  d  =  Tbat.  With  regard  to  the  local  heat  transfer 
in  channel,  the  heat  carried  away  by  airflow  is  equivalent  to  the 
heat  conducted  through  TBL.  Thus,  the  local  heat  transfer  coeffi¬ 
cient  h  in  the  thermal  entrance  region  can  be  expressed  as 


1,  =  -^ 


ey 


y= o 


/Otat-Talr)  =  -AairC/(rbat-rair)  =  ^ 


(12) 


h  keeps  constant  in  the  thermally  developed  region,  and  can  be 
determined  by  substituting  the  local  thickness  of  TBL  at  the  inter¬ 
face  of  the  thermal  entrance  and  thermally  developed  region  into 
Eq.  (12).  Distributions  of  the  local  thickness  of  TBL  <5(n,z )  and  local 
temperature  of  main  airflow  Tair(n,z)  in  the  cooling  channel  n  can  be 
expressed  as  follows: 


<$(n,z)  =  4.518Rez“1/2Pr“1/3z 


(13) 


Tair(n,z)  =T 


amb  ' 


nch 


^airPairQair 


J  h(n)[Tbat(n) 


-Tbat(n-l)-2Tair(n)]dz 


(14) 


where  Rez  is  the  local  Reynolds  number,  Rez  =  uairz/vair,  uair  is  the 
velocity  of  main  airflow.  Since  the  airflow  in  the  cooling  channel  of 
battery  pack  is  generally  laminar,  uair  is  directly  proportional  to 
the  pressure  drop  over  channel  [20].  qair  is  the  volume  flow  rate  of 
air  moving  through  channel.  According  to  the  assumption  (2) 
mentioned  in  Section  2.1,  the  pressure  drop  over  channel  is 
proportional  to  the  total  pressure  drop  across  battery  pack,  and  qajr 
is  equal  to  1/(N  +  1)  of  the  total  airflow  rate. 

For  the  air-cooled  thermal  management  system  of  lithium  ion 
battery  pack,  the  total  volume  flow  rate,  r,  and  pressure  drop, 
APsys,  of  the  air  flowing  across  pack  are  determined  by  the  oper¬ 
ating  point  of  cooling  fan.  Generally,  the  fan  operating  point,  which 
is  dependent  on  the  rotational  speed  of  fan  Hfan  and  resistance 
coefficient  of  battery  pack  Rsys,  is  obtained  by  determining  the 
intersection  of  fan  characteristic  curve  and  system  resistance  curve. 
Fig.  5  shows  the  fan  operating  points  at  different  values  of  rifan  and 
Rsys.  As  can  be  seen  from  the  figure,  with  increasing  fan  speed,  both 
flow  rate  and  pressure  drop  at  the  fan  operating  point  are 


increased;  with  increasing  resistance  coefficient,  the  pressure  drop 
at  the  fan  operating  point  is  also  increased  but  the  flow  rate  would 
be  decreased.  Therefore,  the  critical  parameters  affecting  heat 
dissipation  are  the  rotational  speed  of  cooling  fan  and  the  resis¬ 
tance  coefficient  of  battery  pack.  As  the  fan  speed  is  proportional  to 
the  fan  input  power  Pfan,  the  mathematical  relationship  between 
rifan  and  Pfan  can  be  easily  determined  by  fitting  experimental  data. 
On  the  other  hand,  the  resistance  coefficient  is  related  to  all  the 
structural  parameters  that  can  affect  the  pressure  loss  in  battery 
pack,  such  as  the  width,  length  and  flow  area  of  channel,  the 
inclination  angles  of  upper  and  lower  deflector  plate  and  so  on. 
Considering  that  the  main  purpose  of  this  paper  is  to  establish 
a  model  for  predicting  the  heat  generation  and  dissipation  of 
lithium  ion  battery  within  battery  pack,  the  details  of  modeling  the 
nonlinear  and  interactive  effects  of  the  structural  parameters  on 
resistance  coefficient  would  not  be  discussed  here. 

In  this  paper,  the  values  of  Qnir  and  A Psys  at  different  fan  oper¬ 
ating  points  are  experimentally  determined  and  used  as  the  model 
input.  Substituting  the  measured  uajr  and  Q^r  into  Eqs.  (13)  and  (14) 
to  calculate  the  distributions  of  8  and  Tair,  then  the  heat  dissipation 
rate  of  cell  n  can  be  written  as 

t-ch 

Qdis(n)  =Ach  /  {/i(n)[Tbat(n)  -  Tair(n)] 

0 

+  h(n  +  l)[Tbat(n)-Tair(n  +  l)]}dz 
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Fig.  6.  Flow  chart  for  solving  the  electrochemical  variables. 


3.  Solution  procedure 

To  solve  the  distributions  of  lithium  ion  concentration  and 
electrical  potential  in  electrodes,  we  divide  the  domain  shown  in 
Fig.  2  into  M  control  volumes.  Using  Ax  to  denote  the  length  of 
control  volume,  then  the  control  volume  within  which  reaction  can 
occur  is  SeAx.  It  should  be  noted  that,  the  values  of  Ax  for  the 
electrodes  and  separator  should  be  properly  determined  to  ensure 
that  the  electrode/separator  interface  can  be  matched,  and  that  the 
spatial  variation  of  reaction  current  density  in  each  control  volume 
is  small  enough  to  be  neglected. 

Fig.  6  gives  a  flow  chart  to  illustrate  the  procedure  for  solving 
the  species  concentration  and  electrical  potential  in  negative 
electrode.  As  shown  in  the  figure,  the  reasonable  values  of  the 
variables  j,  ce,  cse,  <ps  and  <pe  can  be  obtained  by  iteratively  modi¬ 
fying  the  initial  value  of  electrolyte  phase  potential  X0,  until  the 
relative  error  of  current  density  reaches  the  preset  threshold  value 


Table  1 

Model  parameters  for  the  144  V/6  Ah  lithium  ion  battery  pack. 


Anode 

Separator 

Cathode 

Design  parameters 

<5n>  <5sep>  (cm) 

50  x  10"4 

25.4  x  10"4 

36.4  x  10"4 

ce  (mol  cm-3) 

1.2  x  10"3 

1.2  x  10"3 

1.2  x  10"3 

Csmax  (mol  Cm"3) 

16.1  x  10“3 

— 

23.9  x  10"3 

Rs  (cm) 

1  x  10"4 

— 

1  x  10"4 

se  (cm2  cm-3) 

1.74  x  104 

— 

1.5  x  104 

Kc(fl) 

1.6  x  10"3[E] 

Se  (cm2) 

1.0452  x  104 

Kinetic,  thermodynamic  and  transport  parameters 

«a,  0LC 

0.5 

- 

0.5 

Oo 

0.11[E] 

- 

0.914[e] 

0i 

0.66[e] 

- 

0.42[e] 

Df  (cm2  s"1) 

4.97  x  10-7 

9.19  x  lO"7 

4.93  x  lO"7 

Ds  (cm2  s"1) 

2.0  x  10"12 

- 

3.7  x  lO"12 

k  (cm2-5  mol-0-5  s"1) 

1.764  x  10-6[11] 

- 

3.626  x  10-6[11] 

/ceff  (Scm-1) 

1.21  x  10"3 

2.24  x  10"3 

1.2  x  10"3 

cjeff  (Scm-1) 

0.58 

- 

0.05 

k\ f  (Acirr1) 

Kf  =  2R7-batKeff 

(t“ -i/f) 

0.363 

Physical  parameters  of  battery  pack 

Cbat  (J  kg-1  K1) 

2460[m] 

mbat  (kg) 

0.696[m] 

Ach  (cm2) 

311.75[m] 

Lch  (cm) 

14.5[m] 

Rcon  (KW-1) 

1.711 [E] 

- 

- 

^Manufacturer. 
[E]Estimated. 
[11]Ref.  [11]. 


6q.  m\  in  the  figure  is  the  number  of  control  volumes  within 
negative  electrode  and  is  equal  to  <Sn/Ax.  Variables  in  the  positive 
electrode  and  separator  can  be  solved  in  the  same  way,  but  the 
values  of  j  and  cps  are  set  at  zero  in  the  range  of  5n  <  x  <  L  -  <5P  as 
there  is  no  reaction  occurred  in  the  separator.  With  the  calculation 
results  of  j,  ce,  <ps  and  <pe,  the  total  heat  generation  rate,  Qgen,  of  cell 
n  at  time  t  can  be  calculated  by  Eqs.  (2),  (7)  and  (8), 
Qgen(Tft)  =  Qr  +  Qo  +  Qp- 

The  distributions  of  TBL  thickness,  heat  transfer  coefficient  and 
air  temperature  in  cooling  channel  can  be  obtained  by  following 
steps:  (a)  substitute  the  air  velocity  ua jr  measured  at  fan  operating 
point  into  Eq.  (13)  to  calculate  the  local  TBL  thickness  <5;  (b)  utilize 
5  to  calculate  the  local  heat  transfer  coefficient  h  by  Eq.  (12);  (c) 
substitute  h  and  the  airflow  rate  Qair  measured  at  fan  operating 
point  into  Eq.  (14)  to  compute  the  local  air  temperature  Tair.  Based 


temperature  sensors 


Fig.  7.  Lithium  ion  battery  pack  and  surface  mounted  thermocouples. 
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Table  2 

Operating  conditions  of  lithium  ion  battery  pack. 


Pfan  (W) 

Pdis  (kW) 

Pcha  (kW) 

Pamb  (°C) 

Case  1 

0 

5 

3 

29.5 

Case  2 

0 

10 

6 

31.3 

Case  3 

1.33 

10 

6 

30.2 

Case  4 

2.43 

10 

6 

27.9 

on  the  calculated  h  and  Tair*  the  cell  heat  dissipation  rate  QdiS(t,n) 
can  be  determined  by  Eq.  (15).  And  then,  combining  Qgen(t,n)  and 
Qdis(t.n),  the  temperature  change  of  cell  n  can  be  deduced  from 
Eq.(l). 

4.  Results  and  discussion 

The  144  V/6  Ah  lithium  ion  battery  pack  given  in  Fig.  7  is  used  as 
the  computational  domain.  As  shown  in  the  figure,  thermocouple 
temperature  sensors  are  mounted  on  the  cell  surfaces  to  measure 
their  temperature  in  real  time.  The  values  of  model  parameters  for 
the  battery  pack  are  given  in  Table  1.  In  the  table,  the  parameter 
values  without  superscript  are  determined  from  Ref.  [7],  and  the 
others  are  determined  from  Ref.  [11],  battery  manufacturer  or 
estimated  by  fitting  the  model  simulation  results  to  experimental 
data.  For  the  modeling  medium,  the  height  and  length  of  the 
control  volume  are  12  cm  and  2  x  10  4  cm  respectively,  and  the 
time  step  size  used  in  the  simulation  is  0.1  s. 

4.1.  Model  validation 

In  order  to  fully  validate  the  effectiveness  of  the  proposed 
thermal  model,  the  predicted  and  measured  cell  temperatures  are 
compared  under  four  different  operating  conditions  given  in 
Table  2.  All  the  operating  conditions  are  composed  by  a  series  of 
charge-discharge  cycles.  Every  cycle  incorporates  a  10  s  constant 
power  discharge,  a  40  s  rest  and  a  10  s  constant  power  charge, 
which  is  similar  to  the  HPPC  test  profile  specified  in  the  battery  test 
manual  [21  ].  PdiS,  Pcha  and  Pfan  in  the  table  represent  the  discharge 
power  and  charge  power  of  battery  pack,  and  the  input  power  of 
cooling  fan  respectively.  Pfan  =  0  indicates  the  battery  pack  operates 
under  natural  convection  condition  and  the  local  heat  transfer 
coefficient  can  be  fixed  at  5  W  m-2  K-1  along  the  whole  cooling 
channel  [22].  At  the  beginning  of  the  charge-discharge  cycling 
experiments  (when  t  =  0),  the  SOC  of  battery  pack  is  50%  and  all  the 
cell  temperatures  are  equal  to  the  corresponding  ambient 
temperature  shown  in  Table  2. 

Figs.  8  and  9  show  the  comparisons  between  the  predicted 
cell  temperature  Tt>at  and  measured  cell  temperature  ?exp  under 
the  four  operating  conditions.  Considering  that  the  temperature 
distribution  is  almost  symmetrical  with  respect  to  the  center  line 
of  battery  pack,  the  figures  only  give  the  first  10  cell  tempera¬ 
tures.  It  can  be  seen  that  the  simulation  results  agree  well  with 
the  experimental  results,  and  that  the  maximum  deviation 
between  the  predicted  and  measured  temperatures  is  less  than 
1  °C,  which  convincingly  proves  that  the  proposed  thermal 
model  can  accurately  simulate  the  dynamic  changes  of  cell 
temperature  when  lithium  ion  battery  pack  operates  under 
natural  and  forced  convection  conditions.  The  rising  rate  of  cell 
temperatures  can  be  efficiently  reduced  by  increasing  the  input 
power  of  cooling  fan,  but  the  temperature  uniformity  across  the 
battery  pack  is  also  deteriorated.  Therefore,  the  thermal 
management  system  of  lithium  ion  battery  pack  should  be 
properly  designed  to  balance  the  temperature  rise  and  unifor¬ 
mity  of  temperature  distribution. 
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Fig.  8.  Comparisons  of  Tex p  with  rbat  under  natural  convection  condition. 

4.2.  Effects  of  current  rate  and  SOC  on  heat  generation  rate 


Fig.  10  shows  the  variations  of  heat  generation  rate  when  the 
battery  pack  with  SOC  of  50%  is  charged  and  discharged  at  different 
current  rate.  Qrev  and  Qirr  in  the  figure  represent  the  predicted 
generation  rates  of  reversible  and  irreversible  heat  of  a  single  cell 
respectively,  Qrev  =  Or,  Qirr  =  Qo  +  Qp-  As  expressed  by  Eq.  (2),  the 
generation  rate  of  reaction  heat  is  proportional  to  the  current 
density,  while  according  to  the  Joule’s  law  the  generation  rate  of 
ohmic  heat  is  proportional  to  the  square  of  current  rate.  Therefore, 
as  shown  in  the  figure,  Qrev  increases  at  an  almost  constant  rate 
with  increasing  current  rate,  but  the  increasing  rate  of  Qirr  is 
increased  gradually. 

Fig.  11  presents  the  variations  of  Qrev  and  Qirr  when  the  battery 
pack  with  different  SOC  is  charged  and  discharged  at  10  C  current 
rate.  It  shows  that  the  SOC  change  affects  QreV  obviously  but  has  no 
significant  influence  on  the  Qirr.  Since  the  0U/07bat  is  strongly 
dependent  on  SOC  as  shown  in  Fig.  3,  the  reaction  heat  generation 
rate  varies  significantly  with  SOC.  Due  to  the  exchange  current 
density  is  closely  related  to  SOC  as  expressed  by  Eq.  (10),  the 
polarization  heat  generation  rate  also  changes  with  SOC  for  a  given 
current  rate.  However,  the  SOC  change  has  no  effect  on  the 
generation  rate  of  the  ohmic  heat  contributed  by  the  contact 
resistance  between  electrode  and  current  collector,  which  accounts 
for  the  biggest  ratio  of  the  total  irreversible  heat.  Thus,  the  change 
of  reaction  heat  with  SOC  is  well  reflected  in  the  variation  curve  of 
Qrev,  while  the  change  of  polarization  heat  generation  rate  with  SOC 
is  not  obvious  in  the  variation  curve  of  Qjrr. 
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Fig.  9.  Comparisons  of  Tex p  with  rbat  under  forced  convection  condition. 
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Besides,  as  can  be  seen  from  the  Fig.  10  and  Fig.  11,  the  variation 
curve  of  Qirr  during  charge  is  almost  completely  coincided  with  the 
one  during  discharge,  and  the  variation  curves  of  Qrev  are 
symmetric.  When  the  battery  pack  with  the  same  SOC  is  charged 
and  discharged  at  the  same  current  rate,  the  number  of  electrons 
and  lithium  ions  collided  with  transport  media,  and  the  change 
extents  of  the  activation  energy  barriers  needed  to  be  overcome  for 
intercalation  and  deintercalation  reaction  are  almost  the  same,  but 
the  entropy  changes  of  electrodes  during  charge  is  exactly  opposite 
to  the  one  during  discharge.  Hence,  the  generation  rates  of 
reversible  and  irreversible  heat  during  charge  are  almost  equal  to 
that  during  discharge,  but  the  signs  of  the  reversible  heat  genera¬ 
tion  rates  are  opposite.  Opposite  sign  means  the  effect  of  reversible 
heat  on  cell  temperature  is  different,  i.e.  if  the  reversible  heat 
generated  during  charge  raised  the  cell  temperature,  then  that 
generated  during  discharge  would  reduce  the  cell  temperature. 

4.3.  Effect  of  resistance  coefficient  on  heat  dissipation  rate 

According  to  the  discussion  presented  in  Section  2.3,  the  resis¬ 
tance  coefficient  of  battery  pack  is  varied  with  the  structural 
parameters  of  cooling  channel,  and  varying  resistance  coefficient 
could  change  the  total  airflow  rate  Qair  for  a  given  Pfan.  With  the 
help  of  the  characteristic  curve  of  cooling  fan,  we  can  use  the 


change  of  Qair  to  describe  the  variation  degree  of  the  resistance 
coefficient  of  battery  pack.  For  example,  if  Qairo  denotes  the  total 
airflow  rate  of  the  144  V/6  Ah  battery  pack  with  its  original  resis¬ 
tance  coefficient,  then  Qair  <  Qairo  means  the  resistance  coefficient 
of  the  battery  pack  is  increased,  and  the  smaller  the  Qair/Qairo,  the 
higher  the  resistance  coefficient. 

Utilizing  given  Qair  and  Pfan  to  determine  the  fan  operating  point, 
we  can  calculate  the  heat  dissipation  rate  for  the  battery  pack  with 
varied  resistance  coefficient.  Fig.  12  presents  the  variations  of  heat 
dissipation  rate  of  a  single  cell  with  resistance  coefficient  at 
different  values  of  Pfan.  The  cell  temperature  and  its  ambient 
temperature  are  equal  to  50  and  25  °C  respectively.  As  shown  in  the 
figure,  there  always  exists  an  optimal  resistance  coefficient  for  the 
battery  pack  such  that  the  heat  dissipation  rate  can  reach 
a  maximum  value.  This  can  be  explained  by  that,  for  a  given  Pfan, 
battery  pack  with  a  lower  resistance  coefficient  corresponds  to 
a  smaller  air  velocity  in  cooling  channel,  which  can  result  in  a  lower 
heat  transfer  coefficient  h  as  expressed  by  Eqs.  (12)  and  (13),  while 
battery  pack  with  a  higher  resistance  coefficient  corresponds  to 
a  smaller  airflow  rate  in  the  channel,  which  would  lead  to  a  higher 
air  temperature  Tair  as  presented  in  Eq.  (14).  Therefore,  according  to 
Eq.  (15),  too  low  or  too  high  resistance  coefficient  is  not  beneficial 
for  dissipating  heat  from  battery  pack.  It  can  also  be  seen  in  the 
figure,  the  optimal  resistance  coefficient  is  increased  when  the  fan 
input  power  increases.  This  is  because  the  impact  of  decreasing  h  on 
heat  dissipation  is  stronger  than  that  of  increasing  Tair  under  big 
input  power  condition.  Thus  it  can  be  concluded  that  when 
a  cooling  fan  with  big  input  power  is  installed  to  provide  cooling 
airflow,  the  thermal  management  system  of  battery  pack  should  be 
designed  to  have  a  high  resistance  coefficient,  in  order  to  make  full 
use  of  the  cooling  ability  of  the  fan. 

4.4.  Effects  of  reversible  heat  on  model  accuracy 

The  lithium  ion  battery  pack  for  pure  electric  vehicle  application 
is  often  charged  or  discharged  continuously,  while  the  one  for 
hybrid  electric  vehicle  application  is  always  cycled  about  a  rela¬ 
tively  fixed  SOC.  In  order  to  comprehensively  evaluate  the  effects  of 
reversible  heat  on  the  model  accuracy,  the  cell  temperatures  are 
predicted  with  and  without  considering  reversible  heat  as  heat 
source,  and  compared  with  the  measured  cell  temperature  under 
different  operating  conditions. 


Fig.  11.  Variations  of  Qrev  and  Qjrr  with  SOC. 


Fig.  12.  Variations  of  QdiS  with  resistance  coefficient. 
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Fig.  13.  Comparisons  of  rexp  with  Tbati  and  Tbat2  during  continuous  charge  and 
discharge,  (a)  constant  current  charge,  (b)  constant  current  discharge. 

Fig.  13  shows  the  comparisons  of  Texp  with  Tbati  and  Tbat2  during 
constant  current  charge  and  discharge  of  1  C  and  5  C  current  rate. 
Tbati  is  the  predicted  cell  temperature  when  Qgen  =  Qirr  +  Qrev  and 
Tbat2  is  the  predicted  cell  temperature  when  Qgen  =  Qirr-  It  can  be 
seen  that  the  reversible  heat  has  significant  effect  on  cell  temper¬ 
ature  when  the  battery  pack  is  charged  or  discharged  continuously, 
and  the  effect  is  especially  obvious  under  low  current  condition. 
This  can  be  explained  by  the  curves  shown  in  Fig.  10.  As  shown  in 
the  figure,  the  generation  rates  of  reversible  and  irreversible  heat 
are  of  the  same  order  under  low  current  condition,  so  the  effect  of 
reversible  heat  on  cell  temperature  must  be  greater  than  that  under 
high  current  condition,  when  the  irreversible  heat  has  become  the 
dominant  heat  source. 

Fig.  14  depicts  the  comparisons  of  Texp  with  Tbati  and  Tbat2  when 
the  battery  pack  with  initial  SOC  of  50%  operates  under  charge- 
discharge  cycles  of  5  C  and  10  C  current  rate.  As  can  be  seen  from 
the  figure,  the  reversible  heat  has  little  influence  on  cell  tempera¬ 
ture  under  this  condition.  This  is  because  the  reversible  heat 
alternately  heat  and  cool  the  cell  when  the  battery  pack  is  cycled 


about  a  relatively  fixed  SOC,  the  effect  of  reversible  heat  on  cell 
temperature  is  eliminated  to  some  degree. 

5.  Conclusions 

(1)  The  proposed  theoretically  based  thermal  model  is  capable  of 
simulating  the  heat  generation  and  dissipation  of  lithium  ion 
battery  pack.  By  charge-discharge  cycling  experiments  under 
natural  and  forced  convection  conditions,  the  accuracy  of 
predicted  cell  temperatures  is  validated  to  within  1  °C. 

(2)  The  simulation  results  indicate  that,  the  change  of  SOC  has 
a  significant  effect  on  the  reversible  heat  generation  rate  but 
has  almost  no  influence  on  the  irreversible  heat  generation 
rate,  and  the  generation  rates  of  reversible  and  irreversible  heat 
during  charge  are  almost  equal  to  that  during  discharge  with 
the  same  SOC  and  current  rate,  but  the  impact  of  reversible 
heat  on  battery  temperature  is  opposite. 

(3)  The  predictions  also  show  that,  for  enhancing  heat  dissipation 
with  a  given  input  power  of  cooling  fan,  there  always  exists  an 
optimum  value  for  the  resistance  coefficient  of  battery  pack. 
Although  increasing  the  input  power  of  fan  would  raise  the  cost 
of  thermal  management,  the  battery  pack  is  allowed  to  have 
a  high  resistance  coefficient,  which  is  helpful  in  improving  the 
compactness  and  volumetric  power  density  of  battery  pack. 

(4)  The  comparisons  between  the  predicted  and  measured  battery 
temperatures  indicate  that,  the  reversible  heat  has  significant 
influence  on  battery  temperature  during  continuous  charge 
and  discharge,  especially  under  low  current  rate,  but  the 
influence  can’t  be  observed  when  the  battery  pack  is  cycled 
about  a  relatively  fixed  SOC. 

References 

[1]  J.  Robert  Selmana,  Said  Al  Hallaja,  Isamu  Uchidab,  Y.  Hirano,  J.  Power  Sources 
97-98  (2011)  726-732. 

[2]  Chanwoo  Park,  Arun  K.  Jaura,  J.  Engines  112  (2003)  1835-1842. 

[3]  Chris  Mi,  Ben  Li,  Derrick  Buck,  Naoki  Ota,  Vehicle  Power  and  ProPulsion 
Conference,  IEEE,  Texas,  2007,  PP.  107-111. 

[4]  D.  Bernardi,  E.  Pawlikowski,  J.  Newman,  J.  Electrochem.  Soc.  132  (1985)  5-12. 

[5]  S.C.  Chen,  C.C.  Wan,  Y.Y.  Wang,  J.  Power  Sources  140  (2005)  111-124. 

[6]  Zhuqian  Zhang,  Li  Jia,  Nan  Zhao,  Lixin  Yang,  J.  Therm.  Sci.  20  (2011)  570-575. 

[7]  Kandler  Smith,  Chao-Yang  Wang,  J.  Power  Sources  160  (2006)  662-673. 

[8]  Lin  Rao,  John  Newman,  J.  Electrochem.  Soc.  144  (1997)  2697-2704. 

[9]  Xiongwen  Zhang,  Electrochim.  Acta  56  (2011)  1246-1255. 

[10]  Dong  HyuP Jeon,  Seung  Man  Baek,  Energ.  Convers.  Manage.  52  (201 1 )  2973-2981 . 

[11]  Yonghuang  Ye,  Yixiang  Shi,  Ningsheng  Cai,  Jianjun  Lee,  Xiangming  He, 
J.  Power  Sources  199  (2012)  227-238. 

[12]  W.B.  Gu,  C.Y.  Wang,  J.  Electrochem.  Soc.  147  (2000)  2910-2922. 

[13]  ChristoPhe  Forgez,  Dinh  Vinh  Do,  Guy  Friedrich,  Mathieu  Morcrette, 
Charles  Delacourt,  J.  Power  Sources  195  (2010)  2961-2968. 

[14]  Noboru  Sato,  J.  Power  Sources  99  (2001)  70-77. 

[15]  Vilayanur  V.  Viswanathan,  Daiwon  Choi,  Donghai  Wang,  Wu  Xu,  Silas  Towne, 
RalPh  E.  Williford,  Ji-Guang  Zhang,  Jun  Liu,  Zhenguo  Yang,  J.  Power  Sources 
195 (2010) 3720-3729. 

[16]  Zhiyu  Jiang,  Jie  Zhang,  Liangjun  Dong,  Jihua  Zhuang,  J.  Electroanal.  Chem.  469 
(1999)  1-10. 

[17]  Marc  Doyle,  John  Newman,  J.  Electrochem.  Soc.  143  (1996)  1890-1903. 

[18]  Weifeng  Fang,  Ou  Jung  Kwon,  Chao-Yang  Wang,  Int.  J.  Energy  Res.  34  (2010) 
107-115. 

[19]  John  H.  Lienhard  IV,  John  H.  Lienhard  V,  A  Heat  Transfer  Textbook,  third  ed., 
Phlogiston  Press,  Cambridge  Massachusetts,  2006,  PP.  286-288,  318-322. 

[20]  R.W.  Fox,  A.T.  McDonald,  P.J.  Pritchard,  Introduction  to  Fluid  Mechanics,  sixth 
ed.,  John  Wiley  &  Sons,  New  York,  2004,  PP.  337-343. 

[21]  DOE/ID-1 1069FreedomCAR  Battery  Test  Manual  for  Power-Assist  Hybrid 
Electric  Vehicles,  Idaho  National  Engineering  and  Environmental  Laboratory, 
2003,  PP.  89-93. 

[22]  A.A.  Pesaran,  S.  Burch,  M.  Keyser,  Proceeding  of  the  4th  Vehicle  Thermal 
Management  Systems  Conference  and  Exhibition,  London,  UK,  1999. 

Nomenclature 

ACh'.  heat  transfer  surface  area  of  cooling  channel  (m2) 

Cbat:  sPecific  heat  caPacity  (J  kg-1 1<-1) 

ce:  lithium  ion  concentration  in  electrolyte  Phase  (mol  cm-3) 

cs:  lithium  ion  concentration  at  the  radial  Position  r  (mol  cm-3) 


164 


C.  Zhu  et  al.  /  Journal  of  Power  Sources  223  (2013)  155-164 


cse:  lithium  ion  concentration  at  the  electrode/electrolyte  interface  (mol  cm-3) 
Csmax •'  maximum  lithium  ion  concentration  in  solid  particles  (mol  cm-3) 

Ds:  diffusion  coefficient  of  solid  phase  (cm2  s-1) 

F:  Faraday’s  constant,  96487  C  mol-1 
h:  local  heat  transfer  coefficient  (W  m-2  K-1) 
io:  exchange  current  density  (A  cm-2) 

I:  applied  current  (A) 

j:  bulk  density  of  reaction  current  in  electrode  (A  cm-3) 
k:  kinetic  rate  constant  (cm2-5  mol-0-5  s-1) 

LCh:  channel  length  (m) 

mba t:  battery  mass  (kg) 

ne:  number  of  electrons  transferred,  1 

Pfan :  input  power  of  cooling  fan  (W) 

Pr:  Prandtl  number  of  air,  0.7023 

qair:  volume  flow  rate  of  air  moving  through  channel  (m3  s-1) 

Qair-'  total  volume  flow  rate  of  air  across  battery  pack  (m3  s-1) 

Qcon-'  heat  transfer  rate  between  two  adjacent  cells  (W) 

Qdis-'  heat  dissipation  rate  (W) 

Qgen'-  heat  generation  rate  (W) 

Qo-'  ohmic  heat  generation  rate  (W) 

Qp.-  polarization  heat  generation  rate  (W) 

Or.'  reaction  heat  generation  rate  (W) 

r:  radial  coordinate  along  active  material  particle  (cm) 

R:  universal  gas  constant,  8.314  J  mol-1  K-1 

Rc:  contact  resistance  between  current  collector  and  electrode  (Q) 

RCOn'-  thermal  resistance  of  the  material  located  between  two  cells  (K  W-1) 

Rs:  radius  of  the  solid  particles  (cm) 

Rsei:  film  resistance  at  the  electrode/electrolyte  interface  (Q  cm2) 

Rsys:  resistance  coefficient  of  battery  pack 

se:  active  surface  area  per  electrode  unit  volume  (cm2  cm-3) 

Se:  active  surface  area  of  electrode  (cm2) 
t:  operating  time  of  battery  pack  (s) 
t°:  transference  number  of  lithium  ion 
Tair local  temperature  of  main  airflow  (K) 

Tamb:  ambient  temperature  (K) 

That •'  battery  temperature  (K) 

Tiay:  air  temperature  in  thermal  boundary  layer  (K) 


Uair:  velocity  of  main  airflow  (ms-1) 

U:  equilibrium  electrode  potential  (V) 
va/r;  kinematic  viscosity  of  air  (m2  s-1) 

APsys:  total  pressure  drop  of  the  air  flowing  across  battery  pack  (Pa) 
AS:  entropy  change  of  electrode  material  (J  mol-1  I<-1) 

Greek  letters 

aa,  olc:  anodic  and  cathodic  transfer  coefficients 
5:  local  thickness  of  thermal  boundary  layer  (m) 

5n:  thickness  of  negative  electrode  (cm) 

5P:  thickness  of  positive  electrode  (cm) 

5sep:  separator  thickness  (cm) 
ee:  electrode  porosity 
p:  activation  overpotential  (V) 

90:  stoichiometry  at  0%  SOC 
6y.  stoichiometry  at  100%  SOC 

Keff:  effective  ionic  conductivity  of  electrolyte  phase  (S  cm-1) 

:  effective  ionic  diffusional  conductivity  of  electrolyte  phase  (A  cm 
Aair:  thermal  conductivity  of  air  (W  m-1  I<-1) 
p-e/jf.  effective  electronic  conductivity  of  solid  phase  (S  cm-1) 
cpe:  electrical  potential  of  electrolyte  phase  (V) 

(ps:  electrical  potential  of  solid  phase  (V) 

Subscripts  and  superscripts 

air:  cooling  air 

bat:  lithium  ion  battery 

ch:  channel 

dis:  dissipation 

e:  electrolyte  phase 

eff:  effective 

fan:  cooling  fan 

gen:  generation 

n:  negative  electrode 

p:  positive  electrode 

s:  solid  phase 

sep:  separator 

sys:  battery  pack  system 


